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AN ITERATIVE TRANSFORMATION PROCEDURE FOR NUMERICAL SOLUTION OF FLUTTER 
AND SIMILAR CHARACTERISTIC-VALUE PROBLEMS ! 


By Myron L. GOSSARD 


SUMMARY 


in tteratice transformation procedure suggested by H. 
Wielandt for numerical solution of flutter and similar charac- 
teristic-calue problems is presented. Application of this pro- 
ceture to ordinary natural-vibrativn problems and to flutter 
problems is shown by numerical examples. Comparisons of 
computed results unih experimental values and with results 
obtuined by other methods of analysis are made. 


INTRODUCTION 


Existing methods of flutter analysis include the 
representative-section method, generalized-coordinate meth- 
uds, matrix methods, and operational methods. The present 
report presents an iteration procedure for analysis of 
flutter and simuar characteristic-value problems. 

For ordinary natural-vibration problems, iterative pro- 
cedures of the Stodola type (references 1 and 2) are suitable 
for finding the fundamental and higher-order natural modes 
and frequencies. The higher-order solutions are obtained 
through use of the orthogonality relations that exist among 
the natural modes. 

Orthogonality relations analogous to those that exist in 
ordinary vibration problems can be found in the flutter 
problem only by introduction of the so-called “adjoint” 
problem. (This additional step is unnecessary in ordinary 
vibration problems by virtue of the fact that they are 
"self-adjoint.") Wielandt has suggested an iterative trans- 
formation procedure (reference 3) which is well-suited to 
the flutter problem in that it avoids the need of orthogo- 
nality relations and hence does not require consideration of 
the adjoint problem. The iterative transformation pro- 
cudure can also be applied to ordinary natural-vibration 
problems with less labor than is generally required in the 
extended Stodola procedure. 

Because both the original and translated works of 
Wielandt are difficult to follow, an explanation of the idea 
of the iterative transformation procedure is given in the 
present report and the application of the procedure to 
ordinary natural-vibration problems and to flexure-torsion 
flutter problems is shown in numerical examples. Com- 
parisons of computed results with experimental values and 
with results obtained by other methods of analysis are 
also made. 
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SYMBOLS 


flexural stiffness 

torsional stiffness 

spanwise coordinate with origin at root of wing 

complex representation of amplitudes and phases 
of translation of elastic axis in harmonic motion 

complex representation of amplitudes and phases 
of rotation about elastic axis in harmonie motion 

coupled mode (y,$) 

complex coeffieients of y wbich, when multiplied 
by y, give complex representation of amplitudes 
and phases of aerodynamic and inertia forces 
associated with translational component of 
barmonic motion 

complex coefficients of 6 which, when multiplied 
by $, give complex representation of amplitudes 
and phases of aerodynamic and inertia forces 
associated with rotational component of har- 
monie motion 

complex coefficients of y which, when multiplied 
by y, give complex representation of amplitudes 
and phases of aerodynamic and inertia torques 
associated with translational component of har- 
monie motion 

complex coefficients of $ which, when multiplied 
by é, give complex representation of amplitudes 
and phases of aerodynamic and inertia torques 
associated with rotational component of har- 
monic motion 

structural-damping coefficients associated with 
flexure and torsion, respectively (see appendix B) 

coefficient of artificial damping (may be either 
positive or negative) 

reduced frequency (bufo) 

frequency of harmonic motion 


1 — 


characteristic value ( i 





length of semichord of wing 

length of cantilever wing from root to tip 

mass ratio (y/mpb?) Î 
velocity of air relative to wing 

distributed mass of wing per unit length of span 
mass density of air 


! Supersedes NACA TN 2316, “An Iterative Transformation Procedure for Numerical Solutioa of Flutter and Similar Characteristle-Value Problems" by Myron L. Gossard, 1951. 
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a distance between midchord axis and elastic axis 
in terms of local semichord, positive when 
clastic axis is behind midchord axis 

“ distance between elastic axis and gravity axis of dis- 
tributed mass of wing in terms of local semichord, 
positive when gravity axis is behind elastic axis 

r radius of gyration of distributed mass of wing 
about elastic axis in terms of local semichord 

G transcendental functions of k (see reference 4) 
time 


F, 

t 

Tu eigenvalue factor ( Cr — 1 ) 
R : 

A 


ratio of complex constants 
length; in numerical solutions, distance between 
specific adjacent stations of wing 


p applied force 

q applied torque 

V shear 

Af bending moment 

et curvature 

8 slope of elastic axis 

T twisting moment 

6 angle of twist 

Subscripts: 

100: true modes or cigenvalues 
a2.a3,a4, .. . transformed modes 

b intermediate derived mode 
HC uus stations | 

R real 

I imaginary 

0 reference value 


51,532,083, . . . 
Superscripts: 
DB... cycles of iteration 

A bar over a symbol indicates a concentrated quantity 
instead of a distributed quantity. 

A prime is used to denote division by «?, 


sweeping functions 


ITERATIVE TRANSFORMATION METHOD OF SOLUTION 
GENERAL FEATURES OF METHOD 


The principle of the iterative transformation procedure is 
similar in form to that of the standard iteration procedure 
for solving characteristic-value problems. Both procedures 
require the determination of the solutions in the order of the 
magnitudes of the eigenvalues, beginning with the funda- 
mental. Both procedures require assumptions of modes, 
integrations which generally must be done numerically, and 
sweeping operations for higher-order-mode determinations. 
The distinguishing features of the iterative transformation 
procedure occur in the determination of solutions higher than 
the fundamental and are as follows: (1) The immediate aim 
is to determine not the true nth mode, as in the standard 
iteration procedure, but a particular linear combination 
composed of all modes from the fundamental to the nth. 
This linear combination is referred to as the transformed nth 
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mode. The transformed nth mode can be made to have 
nodal (zero) points at specified stations of the wing; such a 
feature is highly desirable in numerical work. (2) The 
sweeping operations, which consist of subtractions of lower- 
order-mode shapes from the function obtained by integrating 
the assumed mode, do not employ the orthogonality relations 
as in the standard iteration method but make use of forcing 
functions that, in numerical work, greatly simplify the 
sweeping operations and increase the over-all accuracy of the 
results by making the sweeping operations more consistent 
with the rest of the process, (3) Although the true nth 
eigenvalue is determined directly in the iterative transforma- 
tion procedure, the true nth mode must be computed from 
quantities within the iteration cycle after the transformed 
nth mode is found. 


OUTLINE OF STEPS IN THE PROCEDURE 


The equation of equilibrium of a cantilever beam vibrating 
harmonically in pure flexure is 


d 
de El qv Ê 


or, after integration with proper attention to boundary 


conditions, 
r fz L L "y (dz) o 
y= | | zr |. il yo (dz) (2) 


The solutions of this integral equation are the true natural 
modes (eigenfunctions) yi, Ya, ys, . - . and the corresponding 
natural frequencies (eigenvalues) «ey, os, ws, .... For 
convenience in subsequent discussion, the true modes are 
assumed to be normalized to unity at some position (station 
A) along the beam. | 

The first mode and frequency are assumed to have been 
previously determined by the Stodola process. The iterative 
transformation procedure becomes applicable in the deter- 
mination of the second mode and frequency. As mentioned 
previously, the immediate aim in the iteratiou for the second 
mode is the determination of a linear combination of first 
and second modes which is called the transformed second 
mode. The linear combination y.—y, which has zero 
ordinate at stution À is chosen and defined as the transformed 
second mode to be determined. The iteration for determina- 
tion of this transformed second mode may be described as 
follows: 


(1) A plausible shape ya: for the transformed second mode 
is assumed. This shape must have zero ordinato at station 
A and should satisfy the boundary conditions as closely as 
possible. 


(2) The displacement 


| r fr L L 


resulting from the inertia load yes, (D corresponding to the 
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assumed shape Ya” vibrating harmonically at frequency o 
is calculated. This calculation must usually be done 
numerically with the square of the frequency wz being 
carried along as an undetermined factor. 

(3) A first-mode shape (previously determined) is sub- 
tracted (swept out) from the calculated displacement y, in 
an amount such that the resulting displacement is zero at 
station À. Thus the resulting displacement is 


Jas  — y» — (2 Yi 


(4) The resulting displacement ye? is compared with the 
assumed displacement Y". When the computations are 
numerical, the ratios y,42/y,, are compared at all the sta- 
tions. If the assumed displacement is exactly equal to the 
transformed second mode, the ratios are equal to each other. 
These ratios contain the single unknown w, and the second 
frequency is that value of wa which makes the ratios unity. 

(5) lf the ratios Ya” fat from the first cycle of iteration 
outlined in the four preceding steps are not reasonably the 
same at all stations, the process must be repeated until the 
ratios become reasonably the same. Each new cycle starts 
with the resultant displacement of each preceding cycle. 
The convergence of this process to the second frequency and 
the transformed second mode is proved in appendix À. 

The transformed third mode and the third frequency are 
computed in the following manner. The transformed third 
mode is defined as that combination of the first three natural 
modes which has @ zero ordinate at the same station that was 
used in the transformed second mode (station A) and also a 
zero ordinate at some other station, station B. Thus the 
transformed third mode is defined as 


Ya— U —Y1) 





The iteration is as follows: 

(1) A plausible shape Ya” for the transformed third mode 
is assumed. This shape must have zero ordinates at sta- 
tions 4 and B and should satisfy the boundary conditions as 
closely as possible. 

(2) The displacement 


[f(x 1 L L p 
n= İ f ET [ f Yan Yas (dr) 


is calculated with the square of the frequency e? carried 
along as an undetermined factor. 

(3) The first of two sweeping operations, in which a first- 
mode shape is swept from the displacement y, so as to make 
the resulting displacement at station A zero, is performed. 
This operation gives the displacement 


(on 
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(4) The second sweeping operation, in which a transformed- 
second-mode shape (previously determined) is swept from 
the resulting displacement of step (3) so that the new 
resulting displacement is zero at station B as well as at 
station A, is performed. (This second sweeping operation 
cannot disturb the zero condition at station A established in 
step (3) because the second sweeping function (the trans- 
formed second mode) is identically zero at station A.) Thus, 
the final resulting displacement is 


m Lt | 


(5) Comparisons of the ratios ya? /ya'P at all stations are 
made, and, if they are not reasonably the same, additional 
cycles of iteration are carried out until the ratios become 
reasonably the same. The third frequency is then computed 
from the ratios as explained previously. Convergence of 
this process to the third frequency and the transformed third 
mode is proved in appendix À. 

Frequencies and transformed modes higher than the third 
may be computed by extensions of the process just described. 


PHYSICAL INTERPRETATION OF THE PROCEDURE 


A physical interpretation of the iterative transformation 
procedure can be given. With regard to the transformed 
second mode, for example, the following question may be 
asked: Under what conditions can the beam vibrate in the 
transformed-second-mode shape at the second natural 
frequency? Vibration in shape ya—9s— 3: at frequency o» 
implies an inertia loading yw,*(y.—y,). But if this load is 
substituted in place of yes in the right-hand side of equation 
(2), the result after integration will not be yz—y, but - 


nc İT İsə 6) 


However, if an external (forcing) load of an amount 
y(t — on”) is added to the inertia load, the total load 
y Cog ya— ey) will produce the shape y.—y;. Thus 


| J, E YÜez”ya— a y) (dr) = Yyy (4) 


The inertia and forcing loads are illustrated in figure 1. 
The inertia load acting alone produces a displacement 
(equation (3)) generally different from zero at station A. 
The forcing load produces the displacement 


(25-1) n= IÑ İz gi. [d riada (5) 


This displacement (equal to the sweeping function) has the 
shape of the previously determined first mode and is equal 
and opposite at station À to the displacement due to the 
inertia load; that is, by virtue of the previously assigned 
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normalizations at station A, 


7 2 
E) Qo zn), @ 
Thus the displacement due to the forcing load is completely 
determmed when the displacement due to the inertia load is 
known. The gist of the foregoing analysis is that vibration 
in the transformed-second-mode shape is the response of the 
beam to an oscillatory forcing load of the first-mode shape 
and of frequency equal to the second natural. frequency, 
superimposed on a free vibration of the beam in the second 
natural mode. 
Similar physical interpretations of the iterative transforma- 
tion process for modes higher than the second can be made. 


APPLICATION OF THE PROCEDURE IN ORDINARY COUPLED 
NATURAL-VIBRATION PROBLEMS 


The procedure that has been outlined in a preceding section 
for pure flexure can easily be extended to systems capable of 
simultaneous flexural and torsional displacements. Airplane 
wings belong to the latter class of systems. The only 
distinguishing element in coupled flexural-torsional vibration 
problems is that each natural mode contains, two components, 
the fiexure and the torsion. These components must always 
appear together in a fixed relation to each other. The two 
components must be computed together and must be used 
together. 


(a) 





(a) Transformed second mode: ya” yı yi. 
(b) Inertia load: y«est(j1— pt). 
(c) Foreing load: y(uif—owit) pi. 


FIGURE 1.—lllustration of physical basis of iterative transformation procedure. 


Each coupled mode is a solution of the simultaneous 
differential equations 


d? d? 
FPS El TE yer + bug) (7) 
d 

203 Ly (buy-- Vr) (8) 


Equations (7) and (8), after integration, become (for a 
cantilever beam) 


y= |" f ZI f [ reut bug)(dx)* (9) 


r H 
b il zəl, yaf(buy-L b'r*q)(dz)! (10) 


The solution of the integral equations (9) and (10) for the 
coupled transformed second mode by the iterative trans- 
formation procedure is outlined diagrammatically in figure 2. 
The flexural component. of the displacement for a particular 
step is illustrated in the left-hand side of the figure and the 
torsional component is illustrated at the same level in the 
right-hand side. 


(yit), 


=> s ————————— 


(a) Assumed transformed second mode, 
(b) Intermediate derived mode. 

(c) First-mode sweeping function. 

(d) Derived transformed second mode, 


FIGURE 2.—Illastration of steps in the iterative transformation procedure fur determining 
coupled modes. 


AN ITERATIVE TRANSFORMATION PROCEDURE FOR NUMERICAL SOLUTION OF FLUTTER PROBLEMS 


In the first step, an approximation to a linear combination 
of the true first and second coupled modes is assumed. The 
particular linear combination having zero flexural displace- 
ment at the tip station (station A) is chosen. (For greatest 
numerical accuracy, this nodal point should be chosen in the 
component and at the station where the first coupled mode 
has its maximum numerical value.) The symbols Yya® and 
Pas? are used to designate the flexural and torsional com- 
ponents of this assumed displacement, respectively. In 
general, the magnitude of the torsional component relative to 
the flexural component is difficult to estimate; the most expe- 
dient thing to do is to take one of the components equal to zero. 

The second step is the computation (by numerical inte- 
gration) of the two components of the displacement due 
to the inertia forces ywr (Yat buha) and inertia torques 
qox (buya T D2r202. that are associated with the assumed 
displacement. The result is termed the intermediate derived 
mode. and the symbols x“) and ¢, are used to designate 
its two components. 

The third step is the determination of a sweeping function 
having the shape of the first coupled mode (previously de- 
termined} and a magnitude such that the sum of the inter- 
mediate derived mode and the sweeping function equals zero 
in the flexural component at station A. In algebraic terms, 
the first-mode sweeping function is given by 





a) 
m. ^n 
dnu = Y Pt (12) 


The fourth step is the addition of the intermediate derived 
mode and the first-mode sweeping function to give the derived 
transformed second mode. Thus the two components of the 
derived transformed second mode are 


ye C1 Y (13) 





Ya? — Yp” em 





du = Gy — D 3) $1 (14) 

The ealculation of the ratios Y, 9 /y,,U and ¿22/42 at 
all stations completes the first cycle of iteration. 

Additional cycles are carried out until the ratios at all 
stations in both the flexural and torsional components have 
values that are reasonably thesame. The true second natural 
frequency of the coupled system is then obtained as described 
previously. 

Steady vibration of an airplane wing at zero airspeed is an 
example of coupled natural vibration. The actual numerical 
calculations for the transformed second mode as well as for 
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the first mode and transformed third mode of an airplane 
wing vibrating at zero airspeed are discussed subsequently as 
a special case of flutter. 

The more general equations of airplane flutter at nonzero 
airspeed may be interpreted in such a way that they can 
be solved by a process analogous to that just described for 
coupled natural vibration. 


APPLICATION OF THE ITERATIVE TRANSFORMATION 
METHOD TO FLUTTER 


FLUTTER EQUATIONS 
The differential equations of equilibrium for a wing execut- 
ing simple harmonic motion are 


d? “007 
Ta EI tig) GP + Pè (15) 


—1- 6J (1-499 LB Oy + Qed (16) 


These equations govern & motion represented by 
Y (z, f) — y (utet (17) 
Dr, 0 {ze (18) 


The use of the structural-damping coefficients g, and g, in 
equations (15) and (16) is discussed in appendix B. The ex- 
pressions P,yt+P,y¢@ and Q,y+Q,¢ are the intensities of 
applied force and torque, respectively. For aerodynamic 
and inertia forces and torques due to air flow and distributed 
mass, the P and Q coefficients have values given by the follow- 
ing formulas (rearranged from those in reference 4): 

For P,, 


P,= Pa—tP ry (19) 
in which 
SÊ P BW? /2@ Y NOS 20 
LT (k tita) (Da em 
and 
(5 N2F fr 
Pul C) T (2) at (21) 
For P,, 
P,= Ps, tP (22) 
in which 
bre) 26 2F 2 
P, E [(5-0)5 aret İ(Z ) basê (23) 


and 


THARI os 
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For Qy, 
Q, e Qry—tQry (2 5) 
in which 


sn... 


and 
(NEG) FIOM en 
And for Qg, 
Q. Qro — igre (28) 
in which 


Q -(7) | -Ġ- a° ee a sı a^ HOL o ta? 
(29) 
and 


b NAT 71 1 /1 26 (1 2F 
A Gr) Gi) 1G deu 
(30 


For inertia forces and torques due to concentrated mass, the 
intensities of force and torque are, respectively, 


Pat Ppp lim PtP et (31) 
and 
Qu Q lim Qu Qed (32) 
in which 
E-(x) A) Gs Ga) 
rs 
and 


a= (yt XL) ha (Z) ee (35) 


For a cantilever wing the boundary conditions on the 
displacements are 


(0:-:-(9:-:-(5É) =| Era + 4n) SE] 
uru. (36) 


The differential equations (15) and (16) are now written 
with the eigenvalue & as an explicit factor. Thus equations 
(15) and (16) become 


d? m d £ I F ^ 
El ig) EAP y Py) (87) 
and 


—£ as Kand AQ UE (88) 
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in which the P’ and €” coefficients are equal, respectively, to 
the P and Q coefficients divided by &?. 


FORMULATION OF PSEUDOFLUTTER PROBLEM 


Those solutions a”, (y,¢) of equations (37) and (38) for 
which or is a real and positive (not complex) quantity 
represent the steady harmonic motions of true futter. 
However, because the P and @ coefficients are in general 
complex and because of the presence of structural damping, 
the solutions of equations (37) and (38) will, in general, 
be complex and will include complex eigenvalues a”, As 
in other methods of flutter analysis, the problem is made 
tractable by assuming at the beginning a value of the param- 


eter pote This assumption fixes the values of the P and Q 


coefficients. A real value of & is assumed because » must 
be real and only real values of w can represent flutter. To 
avoid the inconsistence of assumed real values of k and 
obtained complex values of & in the solutions, the problem 
itself is altered by introducing an artificial damping so that 





o? À 
EETA where g; 18 the 
coefficient of artificial damping. Thus the differential 
equations of what may be termed the pseudoflutter problem 


become 


the complex eigenvalue is given by 


d? / 
2, EIO tig T $= D Ti (Py PS $) (39) 


Lai tig) i$ (Q/y+Q'6) (40) 
dz Is da 14 0 FU O 


The value of o” can now be real for any assumed real value 
of k and is therefore the square of the frequency of the steady 
harmonic motion maintained by the artificial-damping 
forces and the naturally present serodynamic, inertia, 
structural, and structural-damping forces. True flutter 
is possible for those special cases in which g, is zero. 

Equations (39) and (40) are similar in form to equations 
(7) and (8) and can be solved by the iterative transformation 
procedure in & way completely analogous to the solution 
of the ordinary problem. The complications introduced 
by the presence of air forces require, however, that a set of 
solutions be obtained for each of several assumed values of £. 
The fact that most of the functions involved are complex 
virtually quadruples the labor as compared with that 
required in the ordinary coupled natural-vibration problem. 


STEPS IN THE ITERATION AS APPLIED TO FLUTTER 


The iteration procedure employs the basic differential 
equations (39) and (40) in their integral forms which, for the 
cantilever wing under consideration, are 


“ol, J. EGF jJ Girt Prada (41 
YO Jo , ETİ tia. 2.5. 


1 = 1 L , f j 
eol, zara, Gr WE) 
- of the 


(42) 


in which O stands for the more convenient form 


eigenvalue. The iteration of equations (41) and (42) 
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follows the same form as the iteration of equations (9) and 
(10). Briefly, the steps are as follows: 

(1) À real value of & is assumed and the values of the 
complex P and Q coefficients are computed. 

(2) An assumption is made for the desired mode y,p. 
(In the first cycle of iteration the assumed mode may be 
real but in the following cycles it will be complex.) 

(3) The complex loadings P,y+2 and Q,y+ Qeg are 
computed. 

(4) The integrations indicated in the equations are 
carried out numerically to get the complex intermediate 
derived mode. 

(3) The sweeping operations are performed by using the 
complex lower-order transformed modes previously deter- 
mined. For convenience in numerical calculations, the 
flexural and torsional components of the complex derived 
(swept) transformed mode are computed in the forms 





4 
Cx BL bo Ky (43) 
and 
1 y L Eal bö” + 
C Ho Eyl, Codo L Kç (44) 


respectively, in which A, and K, are nondimensional com- 
plex functions of the spanwise coordinate z. 

(6) The derived and assumed modes are compared by 
computing their ratios at the stations of the wing. If these 
ratios are not reasonably the same, additional cycles of 
iteration are carried out until the ratios are reasonably 
the same. In the limit (never obtained in practice) the 
ratios will be identical and the proper value of C is that value 
which makes them unity; that is, 


Î Yo L* ” 1 Yo L? Eolo de 

E Ge 5 00 

Û Ho Eolo By _£ Ho Eolo Goda L 57. (45) 
y d 


in which y and ¢ constitute the assumed mode of the limiting 
eycle and the functions in the numerators constitute the 
derived mode of the limiting cycle. 


Equation (45) may be stated in the form 


i . Ya Lt = 
rı Dik) n ËL! (46) 


in which D and H are nondimensional real numbers. Inas- 


much as C is defined as : Tige 





, equation (46) may be written 








Yo L* A Fig 
(H.R t E (47) 
from which the frequency and artificial-damping coefficient 


are obtained as follows: 


Esla Ko 
LA 
a= Y 2 (48) 
H 
=D (4 9) 
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The relative air velocity corresponding to the assumed value 
of Ë is given through the definition of Ë, that is a 


D _ 5 (50) 


NUMERICAL EXAMPLES 


Numerical computations presented in this section illustrate 
the actual application of the iterative transformation pro- 
cedure first to the ordinary natural-vibration problem 
(vibration at zero airspeed) and then to the flutter problem. 
All examples deal with the cantilever wing shown in figure 3. 

The geometric, structural, and mass properties of the 
wing are given in figure 3. A station coordinate system is 
employed for the purposes of the required numerical inte- 
grations. Four stations along the span have been selected 
as indicated in the figure; one of these stations is located at 
the spanwise position of the concentrated mass. The dis- 
tributions of forces and displacements over the span are 
considered to be adequately defined (through interpolation) 
by the forces and displacements at the four selected stations. 
The selection of a system of stations in any problem is im- 
portant because it greatly influences the amount and accu- 
racy of the work to follow. In problems, such as the present 
one, that involve concentrated masses, a station must be 
placed at each concentrated mass because displacements at 
the concentrated masses must be known. (More generally, 
a station must be placed at each discontinuity. Discontinu- 
ities may be present in the distribution of the structural 
stiffness and in the plan form as well as in distribution of 
mass.) The other stations should be equally spaced between 
the discontinuities, and for the system of parabolic inter- 
polation used in the numerical integrations in this report 
there must be a minimum of one station between each 
adjacent pair of discontinuities. The total number of 
stations should be the smallest possible that is consistent 
with the desired accuracy because the calculation effort 
increases rapidly with an increasing number of stations. In 
coupled systems, the number of degrees of freedom allowed 
is twice the number of stations selected; that is the number 
of degrees of freedom in either the flexural or torsional 
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FIGURE 3.—Properties of cantilever wing used İn numerical examples. 
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component of displacement is equal to the number of stations 
employed. Experience has indicated that with parabolic 
approximations results accurate to at least two significant 
figures in the highest mode computed can be obtained by 
employing numbers of stations as follows: For uncoupled 
systems, the number of stations should be two greater than 
the order of the highest mode to be computed; for coupled 
systems, the number of stations should be one greater than 
the order of the highest mode to be computed, with a mini- 
mum of three stations. More than these minimum numbers 
of stations may be required if their use is dictated by suffi- 
ciently many discontinuities. 


ORDINARY COUPLED NATURAL MODES AND FREQUENCIES 


The calculations for the first, second, and third modes at. 
zero airspeed for the wing of figure 3 are shown in tables 1, 
2. and 3, respectively. In this case k= = and the only 
aerodynamic forces are the apparent-mass forces. For 
simplicity, structural damping is disregarded; therefore, all 
quantities entering the problem are real. The numerical 
values of the aerodynamic-inertia force coefficients for 

= co, as well as for other values of k to be used subsequently, 
are given in table 4. 

The first coupled mode is computed in table 1. Table 1 
shows in separate tabulations the flexural and torsional parts 
of the calculation. The first cycle of iteration (part (a) of 
the table) is shown in full detail. Two forms for the tor- 
sional part of the calculations are shown: The. first form 
may be used when the torsional stiffness GJ is constant over 
each bay or over the whole length of the wing; the second 
form, which requires slightly more work, must be used 
when GJ is variable and may be used, as in this case, when 
GJ is constant. The second and third cycles of iteration 
are summarized in parts (b) and (c) of table 1. 

Details of the first cycle of iteration, if the procedure that 
applies only for constant torsional stiffness GJ for the tor- 
sional part of the calculation is used, are as follows: In 
columns 1 of table 1 (a) the two parts y, and 4° of the 
assumed first mode are listed. The torsional component is 
assumed to be zero because it will ultimately be small and is 
dificult to estimate. Columns 2 and 3 are the appropriate 
products of the assumed mode and the distributed-force 
coefficients. Columns 4, which are the sums of columns 2 
and 3, give the two components of the external load which 
correspond to the assumed mode and thearbitrary frequencyw. 
Columns 5 give the concentrated loads (external forces and 
torques) that are equivalent to the distributed loads of 
columns 4. These equivalent concentrations are given in 
columns 5 in terms of the pertinent distances between sta- 
tions ^, and in columns 6 in terms of the reference distance 
X. Formulas used for computing the equivalent concen- 
trations from the distributed loads are given in appendix C. 
Columns 7 and 8 are the appropriate products of the assumed 
mode and the concentrated-force coefficients. Columns 9 
are the total concentrated loads, the sums of columns 6, 7, 
and 8. 
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The flexural and torsional calculations must now be de- 
scribed separately. In column 10 for flexure, the average 
shears in the bays between stations are found by a successive 
summation of the concentrated loads from the tip where 
the shear is zero inboard to the root. In column 11 the incre- 
ments of bending moment are computed by multiplying 
the shears by the bay lengths in terms of X. The bending 
moments of column, 12 are found by a successive summation 
of the increments of bending moment from the tip where the 
bending moment is zero inboard to the root. Column 13. 
gives the distribution of curvature, which is obtained by 
dividing each ordinate of the bending-moment curve by the 
local value of ET (El in this example is constant). Equiva- 
lent concentrated curvatures are now obtained by applying 
to the distributed curvatures the previously used formulas 
for equivalent concentrations. Column 14 gives these 
equivalent concentrations in terms of the distances As, and 
column 15 gives them in terms of the reference distance X. 
The average slopes in the bays are obtained in column 16 
by a successive summation of the concentrated curvatures . 
from the root where the slope is zero outboard to the tip. 
The increments of derived flexural displacement are com- 
puted in column 17 by multiplying the average slopes by the 
bay lengths in terms of Xy. The flexural component y,” of 
the derived mode is obtained in column 18 by a successive 
summation of the increments of displacement from the root 
where the displacement is zero outboard to the tip. Column 
19 gives the ratios at the selected stations of the derived 
flexural component to the assumed flexural component. 

Columns 10 to 15 for torsion are now considered. Column 
10 gives the average twisting moments in the bays of the 
wing and is obtained by a successive summation of the 
concentrated torques of column 9 from the tip where the 
twisting moment is zero inboard to the root. The average 
twists in the bays are computed in column 11 by dividing 
the average twisting moment in each bay by the local value 
of GJ (GJ in this example is constant over the whole span). 
The increments of derived torsional displacement are obtained 
in column 12 by multiplying the average twists by the bay 
lengths in terms of X, The torsional component of the 
derived mode is computed in column 13 by a successive 
summation of the increments of displacement from the root 
where the displacement is zero outboard to the tip. Inas- 
much as the derived displacement of column 13 is in terms 
of GJ, the displacement is converted into terms of EI in 
column 14 so that it may be compared with the assumed 
torsional displacement on the same basis as the assumed 
and derived flexural displacements are compared and so that 
the next cycle may be started with displacement components 
having the same dimensions as the assumed mode of this... 
first cycle. Column 15 normally would contain the ratios 
at the selected stations of the derived torsional component 
to the assumed torsional component, but in the case of table 
1 (a) these ratios are meaningless because the torsional 
component will ultimately be different than was assumed 
in column 1. 
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Before the results of further cycles of iteration for the first 
mode are described, the form that the numerical integration 
for the torsional component must take when GJ is variable 
is described. In the part of table 1 (a) showing the calcula- 
tion for variable GJ, columns 1 to 4 are the same as in the 
calculation for constant GJ. The form of the numerical 
integration changes at column 5. Column 5 consists of 
increments of twisting moment over the bays. "These incre- 
ments are obtained as increments of area beneath the curve 
of distributed torque (column 4). Formulas used for com- 
puting these increments are given in appendix C. In column 
5 the inerements of twisting moment are given in terms 
of the distances X and in column 6 they are given in terms 
of the reference distance Ap. The twisting moments at the 
selected stations due to the distributed torsional loading are 
obtained in column 7 by a successive summation of the incre- 
ments of twisting moment. "The components of external 
concentrated torque are as for constant GJ and are given In 
columns 8 and 9. ‘The applied concentrated torque gives 
twisting moments as shown in column 10. Column 11 is 
the sum of columns 7 and 10 and gives the total twisting 
moments at the selected stations. (Note that in columns 10 
and II there is a discontinuity in twisting moment at the 
station having the mass discontinuity.) Column 12 gives 
the distribution of twist found by dividing column 11 by 
the local value of GJ (GJ being in general not constant). 
The increments of derived torsional displacement are com- 
puted in columns 13 and 14 by applying to the values of 
eolumn 12 the same formulas applied previously to column 
4. The torsional component of the derived mode (columns 
15 and 16) is, except for small computational discrepancies, 
the same as in the previous method, as it should be. 

Two additional cycles of iteration were found to be ade- 
quate for the determination of the first mode and frequency. 
The results of these iterations are shown in parts (b) and (c) 
of table 1. In table 1 (b), for example, columns 1 give the 
two components of the assumed mode of the second cycle, 
which aro obtained by normalizing the derived mode of the 
first cycle to unity in the flexural component at the tip 
station. This normalization is not essential but facilitates 
manipulations and comparisons by keeping the numerical 
values in all cycles within the same range of magnitude. 
Columns 2 give the derived mode obtained by the numerieal 
integration procedure just described. The ratios of derived 
to assumed mode are given in columns 3 for both components 
of displacement. These ratios are seen to be fairly uniform. 
The ratios obtained in the third cycle in table 1(c) are, for 
practical purposes, Identical. The averaging device shown 
in columns 4 of table 1 (c) and below table 1 is adopted 
us a quick and generally quite accurate way of smoothing 
out small discrepancies that remain in the ratios after 
convergence is almost complete. This device, although 
clearly not necessary in the case of table I (c), is useful in 
other cases throughout the numerical examples and is ex- 
plained as follows: The two ratios in columns 4 are obtained 
by considering the flexural and torsional components of the 


displacement separately and then dividing the sum of the 
station values of the derived displacement by the sum of the 
station values of the assumed displacement. When a dis- 
erepancy remains between two ratios of the type in columns 
4, the average of these two is taken as the final value; 
the final value for this case is given in the calculation 
below table 1. This device gives greater weight to the larger 
ordinates and is in that respect similar to other weighting 
procedures such as the energy and least-squares methods 
but is much simpler. If the assumed and derived displace- 
ments contain both positive and negative ordinates, the 
negative ordinates should be changed to positive for the 
purpose of the summations. The remaining calculation 
shown below table 1 gives that value of the arbitrary fre- 
quency w which makes the ratio just computed unity. As 
proved in appendix A, this value of w is the fundamental 
frequency or. 

Table 2 gives the main results of three cycles of iteration 
required to obtain satisfactory approximations of the second 
frequency and the transformed second mode at zero airspeed 
(k= co). Columns 1 of the first cycle (parts (a) of table 2) 
contain the two components Y.” and $4! of the assumed 
transformed second mode. This mode must have one zero 
ordinate (excluding the root ordinates). Although this 
zero ordinate may theoretically be taken at any station, the 
numerical accuracy of the results is greatest if the zero ordi- 
nate is placed at the station and in the component where the 
preceding mode (the first) has its maximum numerical value 
(since the numerical process is such that the larger ordinates 
contain more significant figures than the smaller ordinates). 
Therefore, the zero ordinate of the transformed second mode 
is placed at the tip station in the flexural component, this 
location being designated station A. In the numerical 
values of columns 1, the flexural component ya would 
normally be taken as zero. (The values that are shown are 
estimated from flutter calculations that had previously been 
made for this wing.) Columns 2 give the intermediate 
derived mode obtained by numerical integration. Columns 
3 constitute the first-mode sweeping function. The shape 
of this sweeping function is given by columns 2 of table 1 (c) 
and its magnitude is such as to be equal and opposite to the 
intermediate derived mode at station A. Thus the derived 
transformed second mode (columns 4), which is the sum of 
columns 2 and 3, has zero ordinate in the flexural component 
at station A and a shape comparable to the assumed mode, as 
indicated by the ratios in columns 5. The ratios in the next 
two cycles (parts (b) and (c) of table 2) show marked im- 
provement in uniformity. The final value of the ratio com- 
puted below the table gives, as proved in appendix A, the 
value of the second frequency we, as shown. 


The main results of the iterations to obtain satisfactory 
approximations of the third frequency and the transformed 
third mode at zero airspeed are stated in table 3. Typical 
operations required in & cycle are outlined in table 3 (a). 
Columns 1 give the assumed transformed third mode made 
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up of the two components ya? and ga”. The transformed 
third mode is to have & zero ordinate in the flexural compo- 
nent at the tip station as in the transformed second mode and 


a zero ordinate in the torsional component at the tip station.- 


The location of the second zero ordinate is designated station 
B. To obtain greatest numerical accuracy, the selection of 
the second zero ordinate is governed by the same rule that 
was used for selecting the first zero ordinate, namely, that 
the new zero ordinate should be placed at the station and 
in the component where the preceding mode (the transformed 
second) has its maximum numerical yalue. The numerical 
values thet are shown in columns 1 are estimated from pre- 
vious flutter calculations; the torsional component 440 
would normally be taken as zero. Columns 2 give the inter- 
mediate derived mode, and columns 3 give the first-mode 
sweeping function which, as before, has a magnitude at sta- 
tion A thet is equal and opposite to the intermediate derived 
mode. Columns 4 constitute the transformed-second-mode 
sweeping function which bas a shape given by columns 4 of 
table 2 (c) and a magnitude at station B equal and opposite 
to the sum of the intermediate derived mode and the first- 
mode sweeping function (the sum of columns 2 and 3). The 
derived transformed third mode of the first cycle is the sum 
of columns 2, 3, and 4 and is given in columns 5. ‘The ratios 
in columns 6 are far from uniform. The ratios in the second 
and third cyoles (parts (b) and (c) of table 3) show improve- 
ments in uniformity. The iteration is discontinued at the 
end of the third cycle where tho ratios are about as uniform 
as they can get with the limited number of significant figures 
that are present. The frequency obtained by the smoothing 
device is the third frequency « and has the value shown. 

The patterns leid out in the foregoing examples establish 
the general technique that can be used to obtain zero- 
airspeed modes and frequencies higher than the third. 
Guiding rules for determining the number of selected stations 
to be employed have been given previously. These examples 
also set the basic pattern for the computation of the modes 
and eigenvalues of pseudofiutter and of flutter. 


MODES AND EIGENVALUES OF PSEUDOFLUTTER AND OF FLUTTER 


The operational solution in reference 5 gave for the wing 
under consideration (fig. 3) a reduced_frequency at flutter 
of 0.1443. In order to use this operational solution, this 
same value (k—0.1443) is used in the flutter calculations 
that follow. 

The calculations for the first, second, and third modes at 
k=0.1443 are shown in tables 5, 6, and 7, respectively. 
Aerodynamic-inertia force coefficienta have been computed 
by equations (19) to (35) and their values are given in 
table 4. - Structural damping is disregarded, although a note 
on the method of incorporating structural damping in the 
calculations is made subsequently. 
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Table 5 (a) shows in detail the first cycle of iteration for 
the first mode. The form of the computations is the same 
as that shown previously for the determination of zero- 
airspeed modes. The amount of computation, however, is 
between three and four times that required for zero-airspecd 
modes because of the fact that the functions involved are 
complex and thus must be described by two parts—a 
real part and an imaginary part. Columns 1 and 2 are the 
real and imaginary parts, respectively, of the assumed first 
mode. Asa start, all parts of the assumed mode except the 
real part of the flexural component are taken as zero. Col- 
umns 3 to 6 are the real parts of the products of aerodynamic- 
inertia cocfficients and the assumed mode, and thus their 
sums (columns 7) are the real parta of the distributed load. 
If the expressions for the distributed load are considered, 
this condition is more evident. The distributed forces pro- 
ducing flexure are given by 


(Py —tPr)yettynt+ Pas—tPre)(Oetign=Payet 
Pedes L Pwr Pad il Pyr t Pa ,$r— Py s Pade 
(51) 


The terms of the real part of equation (51) appear in columns 
3 to 6 in the flexural part of table 5 (a); the terms of the 
imaginary part of equation (51) appear in columns 22 to 25 
in the flexural part of table 5 (a). This seperation of real 
and imaginary parts allows the displacement due to each part 
to be computed separately. A similar explanation can be 
made for the quantities in columns 3 to 6 and columns 18 to 
21 in the torsional part of table 5 (a). 

Real and imaginary parts of the concentrated loads that 
are equivalent to the distributed loads are computed as 
explained previously by the formulas of appendix C. These 
values are shown in columns 8, 9, 27, and 28 in the flexural 
part and in columns 8, 9, 23, and 24 in the torsional part. 
The real and imaginary parts of the loads due to the con- 
centrated mass follow next in order, and the total concen- 
trated loads are given in columns 12 and 31 in the flexural 
part and in columns 12 and 27 in the torsional part. The 
average shears, average twisting moments, and bending 
moments are then computed as described previously. 

The remaining parts of the computations in table 5 (a) 
thet are associated. with the real parts of the load are de- 
scribed as follows (the remaining parts that are associated 
with the imaginary parts of the load are similar): Column 
16 in the flexural part gives the distributed curvature duo 
to the real part of the load. This curvature is obtained by 
dividing the ordinates of the real part of the bending-moment 
curve by the local values of the complex flexural stiffness 
EI(1--igj)(1--ig). In these examples, any actual struc- 
tural damping is disregarded: therefore g, is zero. The 
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factor 1-F47,, containing the as yet unknown artificial- 
damping coefficient, combines with & to give the factor 1/C 
.. 
in column 16, € being the arbitrary eigenvalue DCN. 
If the actual structural damping g, is regarded as other than 
zero, the values in column 16 would be computed as follows: 
The real and imaginary parts of the bending moment would 
be combined into the complex bending moment Jfe+iMr. 
This complex bending-moment distribution would then 
be divided by the local values of the complex stiffness to give 


Ma iM; The factor 1+ig, would be carried 


EIU ig) (12-39. 
along in the arbitrary eigenvalue C, and the numerical values 


Mf I+ À, 


of the real part of the quotient ET Fig) would be placed 


in column 16. The imaginary part of the quotient would 
be similarly placed (in column 35) in the calculations asso- 
ciated with the imaginary part of the load. The average 
twists due to the real part of the load are computed in 
column 14 in the torsional part of table 5 (a), and those 
due to the imaginary part of the load must also be computed. 
These calculations follow the same pattern as those just 
explained for the curvatures. The complex torsional stiff- 
ness GJ(1+ig.)(i+ig.) enters in place of the complex 
flexural stiffness. If GJ or gs is variable over a bay length 
or over the whole span, the numerical integration for the 
torsional part of the calculations should be carried out as 
explained in the part of table 1 (a) that deals with variable GJ. 

The numerical integrations are completed in the manner 
already described, and the derived mode is thereby obtained 
in the form of four components of displacement. The 
flexural components are qr? and yii? of columns 21 and 40 
in the flexural part. The torsional components are dış”? and 
du? of columns 17 and 32 in the torsional part. However, 
these components are not actually the real and imaginary 
parts of the flexural and torsional components of the derived 
mode, because each one of them contains the complex factor 
I+ig,. Nevertheless, the complex derived mode is given 
by Yı? üzü” and die Hipu”. 

The complex ratios of the complex derived mode to the 
complex assumed mode are computed in column 41 in the 
flexural part and column 33 in the torsional part. Only 
two of these ratios have actually been computed but they 
are sufficient to indicate the need for further cycles of 
iteration. 

A total of four cycles of iteration (the main results of the 
last three are shown in parts (b), (c), and (d) of table 5) 
was required for satisfactory convergence. In columns 6 of 
tuble 5 (d) and immediately below table 5, the smoothing 
device described previously is applied to obtain the best 
single value of the ratios. The fundamental (first) eigen- 
value is that value of O which makes the ratio unity. Thus 


212183—854——-21 


zy. 


C,—(269.5 —82 21) pr and since C, is defined as 


1 I +19: 
2 aq 4 

Q1 
the frequency and —- damping of the first mode are 
obtained from the real and imaginary parts of the equation 


-— —(269.5—82.2i) ETL E (52) 


The caleulation of these quantities and the corresponding 
airspeed », which is obtained from the relation e are 


shown at the bottom of table 5. 


Tables 6 and 7 show the main results of the iterations to 
obtain the transformed second and third modes for £=0.1443. 
Four eycles of iteration for each mode gave satisfactory 
convergence. The assumed modes of columns 1 and 2 of 
tables 6 (a) and 7 (a) were taken in the forms recommended 
previously in connection with tables 2 and 3. In tables 6 
and 7, the complex intermediate derived modes are given 
by re” Hiyar ™ and der Y 4-146477, the complex first-mode 
sweeping functions, by Yar?” Hiyar™ and due” Flid 
with shapes corresponding to columns 3 and + of table 
5 (d), and the complex transformed-second-mode sweeping 
functions, by ges" Fiya and dur" Hipa with 
shapes corresponding to columns 7 and 8 of table 6 (d). The 
results computed in and below table 7 give for the third 
eigenvalue g,,—0.030 and w;=168.9 radians per second. 
The corresponding sirspeed is v4=390 feet per second. 


COMPUTATION OF TRUE MODES 


Because the critical flutter velocities are given directly by 
the eigenvalues, knowledge of the true modes in flutter 
problems is of no value (at least of no value recognized at 
present). The same statement applies to the transformed 
pseudoflutter modes, with the exception that in the iterative 
method their determination is a necessary adjunct to the 
determination of the eigenvalues. In ordinary problems of 
forced vibration (at zero airspeed), however, the true modes 
are often used with great advantage. For this reason and 
for the sake of completeness of the presentation of the itera- 
tive transformation procedure, the method of determining 
true modes from results of the iterative transformation 
procedure is illustrated in tables 8 and 9. 


The computations in tables 8 and 9 pertain to the same 
wing analyzed in the previous examples. The modes com- 
puted are for k=0.1443. The true third mode as computed 
in table 9 may therefore be compared with the flutter 
mode computed for this wing by the operational method in 
reference 5. 

In table 8, the true second mode is computed as follows 
from functions appearing in the last eycle of iteration for the 
transformed second mode (table 6(d)): Preceding the table 
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proper is the calculation of the eigenvalue factor F wel 
2 


that is needed for computing the true second mode. In the 
terminology of tables 6(d) and 8 and as shown in appendix A, 
the true-second-mode shape is given by 


Yort+tYor=(Yiattyint (Yan? + Yaga") (5 3) 
and 


dan - 1 bar = (dir tt dir) (@nan'? + $4.59) (54) 


in which 


) 
Vuoi o P Ha (55) 
AT! 
and 
dın-İ-14::— dns 177: - (56) 
E 


Columns 1 of table 8 show the key ordinate (Yu 4-195; 9)4 of 
the first-mode sweeping function yy s ® tiyan: ®, doin tidy 
as given in columns 5 and 6 of table 6(d) The key 
ordinate is taken as the largest ordinate (the ordinate at 
station A) for the reason of accuracy cited previously. The 
key ordinate of the first-mode shape yis-- iyi, ret tdi (equal 
to the first terms on the right-hand sides of equations (53) 
and (54)) is shown as the boxed value in column 2 and is 
obtained by dividing the value in column 1 by the eigenvalue 
factor Fi. The other values in columns 2 are obtained by 
using the key ordinate in conjunction with the first-mode 
shape given in columns 3 and 4 of table 5(d). Columns 3 
show the transformed-second-mode shape Zen” Hiyar, 
Par Higor? (equal to the second terms on the right- 
hand sides of equations (53) and (54)) given by columns 
7 and 8 of table 6(d). The sum of columns 2 and 3 which 
is given in columns 4 gives the shape of the true second 
mode Yor+tYor, dartider (equal to the left-hand sides of 
equations (53) and (54)). 

In table 9, the computation of the true third mode 7 


ceeds as follows: The necessary eigenvalue factors F vN 1 


and Fa rii are computed as shown. In the terminology 


of tables 7(d) and 9 and as shown in appendix A, the true- 
third-mode shape is given by 


Van Var (Yue + Yun + (Hier + LY 121) + (Yar + Yan) + 


(#asn' agr) (57) 
and 
dar T 1 dar (duin tipu) + (ier F dion + (boon + idan): c 
(har +t daar) (58) 


in which 


C C 
(2 1) (Yurt im) (E 1 ) (ie VY 122) 
— fo” + yor (59) 


(E = 1) (ue in «(8 = 1) (pier T 16121) 


= bon” +t odur? (60) 
O 
Yar Ye = dı. ua. (61) 
o 1 
Ls (4) 
door tig ——.—.. — o. (69) 
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and Vin t+ Tar, diortidier is to Yaor + Yess; par tipar 8S 
Virt Ur, dirt Pu is to Van -- tar), bar? + ipar” in table 
8. The key ordinates ant? 4-295; 9)4 and (drar + 167%) 
of the first and second sweeping functions appear in columns 1 
and 2 and are taken from columns 5, 6, 7, and 8 of table 
7(d). The key ordinate of the functions Fas(YarH ita), 
Falur tiginn, Which are equal to the second terms on the 
left-hand sides of equations (49) and (60), is computed in 
columns 3 by using the key ordinate of columns 2 in con- 
junction with ordinates at stations À and B in columns 2 
and 3 of table 8 as follows: 


[Fos( yis T L 12D) A] cabreo 
“Kr (Yin ttyina 
(6 


aon ida )a 


The key ordinate of the functions FistVun-b Yu), Fadu sK idu d, 
which are equal to the first terms on the left-hand sides 
of equations (59) and (60), is given in columns 4 and, 
m accordance with equations (59) and (60), is the difference 
between Your? +74”, don tid? of columns 1 and 
Foa(Yierttyier), Fas(Prartibiar) of columns 3. The key ordi- 
nates of the first-mode shapes Yurt+iYur usd idu; and 
ian T Saar, rar tider are shown in columns 6 and 5 and are 
obtained by dividing the values in columns 3 and 4 by the 
appropriate eigenvalue factors. 'The sum of the key ordi- 
nates of columns 5 and 6, shown as the boxed value in col- 
umn 7, is the key ordinate of the total-first-mode shape 
ig Wir, die tidy Which is equal to the sums of the first two 
terms on the right-hand sides of equations (57) and (58). The 
other values in columns 7 are obtained by using the key ordi- 
nate in conjunction with the first-mode shape given in columns 
3 and 4 of table 5(d). The key ordinate of the transformed- 
second-mode shape #aan-+ far, Partida, Which is equal 


aai ua... | [Fast dart tanel table 9 (63) 
table & 
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to the third terms on the right-hand sides of equations 
(57) and (58), is shown as the boxed value in columns 8 and is 
obtained by dividing the value in columns 2 by the eigen- 
value factor F4. The other values in columns 8 are com- 
puted by using the key ordinate in conjunction with the 
transformed-second-mode shape given in columns 7 and 8 of 
table 6(d). Columns 9 show the transformed-third-mode 
shape en” + yar, dar” Hip (equal to the fourth 
terms on the right-hand sides of equations (57) and (58)) 
given bv columns 9 and 10 of table 7(d). The sum of col- 
umns 7, 8. and 9 given in columns 10 gives the shape of 
the true third mode vast tiar, Partida (equal to the left- 
hand sides of equations (57) and (58)). 


TRENDS AND COMPARISONS OF NUMERICAL RESULTS 


Results of the computations shown in the preceding 
section of the report together with results of similar compu- 
tations based on other assumed values of k are given in 
figures 4 to 6. Figures 4 and 5 deal with the wing to which 
the concentrated mass is attached. Figure 6 gives data of & 
similar nature for the same wing without the concentrated 
mass. The computed results obtained by the Rayleigh- 
Ritz and operational methods and the experimental results, 
all of which are given for this wing in references 5 and 6, are 
also recorded in figures 4 to 6. 

In part (a) of figure 4 the solid curves show the variation 
of the artificial-damping coefficient ge with airspeed in cach 
of the first three solutions. For each assumed value of k & 
dashed curve is drawn through points that represent solutions 
for that value of & Part (b) of figure 4 shows in & similar 
way the variation of the frequency w with airspeed and the 
lines of constant values of +. The facts of particular interest 
that are shown by these plots are as follows: 

(1) The true flutter condition is given by the third solution 
for a value of k between 0.1443 and 0.1590 at an airspeed 
almost equal to that found in the experiment. Here the 
computed value of g, is zero. The computed frequency at 
true flutter is also in very close agreement with the experi- 
mental value. 

(2) The operational solution is in good agreement with the 
experimental solution, but the solutions obtained by the 
Rayleigh-Ritz method with three and four modes vary by 
72 percent and 22 percent, respectively, from the solution 
obtained by the operational method. ‘The operational solu- 
tion is theoretically the most exact even though it involves 
summations of finite numbers of terms of infinite series. 
However. as pointed out in reference 5, its use is limited in 
practice to wings of uniform section. In the present example 
the results obtained by the iterative method would be 
expected to be better than the results obtained by the 


Rayleigh-Ritz method because the eight degrees of freedom 
used in the iterative method are much less restrictive than 
the three or four used in the Rayleigh-Ritz method. 
Although exact agreement of the results of any of the com- 
putational methods with the experimental results is not to 
be expected, the better agreement of the iterative solution 
as compared with the operational solution is at first sur- 
prising. On further observation, however, this agreement 
must be credited to a fortunate disposition of the errors 
involved in the iterative method because, in the case of 
figure 6, the relative order of agreement of the operational 
and iterative results with the experimental result is opposite 
to that in figure 4. 

(3) The trends of the solid curves representing the first 
and second solutions in figures 4 (a) indicate that both may 
cross the zero artificial-damping axis at very large airspeeds. 


Experiment! (reference 5) 
Operational solution (reference 5) 
Rayteigh-Ritz: 9 modes frefererce 6) 
Rayleini-Fitz 4 modes (reference 6) 
/feraHarr 4 stations 
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(a) Variation ofartificial damping with airspeed. 
(b) Variation of frequency with airspeed. 


FIGURE 4.—Varistion of artificial damping and frequency with airspeed In first three solutions. 
Wing with concentrated mass, 
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But this conjecture is of no practical interest so long as a 
eurve (the third solution) that crosses at a lower airspeed 
exists. However, the question of whether the curve for 
some solution higher than the third could cross the zero 
artificial-damping axis at an airspeed lower than that at 
which the third solution crosses demands an answer. 

(4) Reasonable assurance that, among all possible solu- 
tions, the curve of third solutions in figure 4 (a) crosses the 
zero artificialdamping axis at a lower airspeed than any 
other is provided by the trends of the curves for constant 
values of k in parts (a) and (b) of figure 4. The curves of k 
sbow that the curve representing the fourth solution will 
most assuredly lie above and to the mght of the solid curves 
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in figure 4 (b) and probably below and to the right of the 
solid curve for the third solution in figure 4 (a). The curves 


of k in figure 4 (b) are straight lines by definition (a =>) 


Prediction of the courses of the curves of k in figure £ (a) 
cannot be made with much certainty. They have a strong 
tendency to proceed to the right, but it is easy to believe that 
upward or downward changes in their directions could take 
place. The curve for the fourth solution, however, would 
probably cross the zero demping axis at a value of r between 
500 and 600 feet per second in figure 4 (a). | 

Figure 5 shows and compares the amplitude and phase 
distributions of modes computed by the iterative transfor- 
mation procedure and by the operational method for the 
wing with a concentrated mass. The first and second modes 
as well as the more important third mode from the iterative 
solution for k=0.1443 are plotted, and . third mode İrom 
the iterative solution for k=0.1590 is also plotted. The 
third modes from the iterative solutions for the two values 
of k agree very well in shape with the flutter mode obtamed 
in reference 5 by the operational method, and the operational 
mode lies between the two iterative modes. Thus the 
agreement of the iterative and the operational methods is 
apain evidenced. | 

Figure 6 is & plot similar to figure 4 but relates to the 
behavior of the wing analyzed in figure 4 if the concentrated 
mass is not present. There is very little similarity in the 
data of the two figures. The most notabie difference is that 
in figure 6 the true flutter mode appears in the second solution 
instead of the third as before and that the flutter speed is 
lower than before. Of interest is the occurrence of almost 
equal eigenvalues in the second and third solutions for 
k=0.50. The flutter speeds given in figure 6 by &ll methods 
of solution, including the Rayleigh-Ritz method, are seen to 
be in substantial agreement. 


CONCLUDING REMARKS 


The report has described the iterative transformation 
method suggested by H. Wielandt and has demonstrated 
the use of the method in an orderly computation of critical 
flutter speeds. Numerical comparisons with solutions ob- 
tained by other methods and with experimental values have 
been made. The applications made in this report show 
promise for future practical use of the method. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY CoMMITTEE FOR ABRONAUTICS, 
LANGLEY FIELD, VA., January 17, 1951. 


APPENDIX A 
ON THE CONVERGENCE OF THE ITERATIVE TRANSFORMATION PROCEDURE 


INTRODUCTION 


The extensive existing literature on the eigenvalue prob- 
lems is concerned almost exclusively with the class known as 
self-adjoint problems, in which the eigenfunctions and eigen- 
values are real. In recent years, non-self-adjoint eigenvalue 
problems have received increasing attention. This class 
includes the flutter problem in which the eigenfunctions and 
elgenvalues are generally complex. The literature referred 
to by Wielandt in reference 3 reveals that the non-self-adjoint 
eigenvalue problem and the transformation method for its 
solution have been given some attention since at least 1928. 
Wielandt's own work constitutes probably the most extensive 
contribution on the subject. 

The discussion on convergence given herein is not con- 
tamed in Wielandt’s work and may be considered a rigorous 
proof if the following assumption is valid: that the equations 
(equations (41) and (42)) for the system (the wing) under 
consideration have an infinite number of solutions that form 
a complete set for any value of the reduced frequency k. 
In the subsequent demonstrations, the validity of expanding 
arbitrary displacement functions in infinite series of eigen- 
functions depends upon the validity of the assumption. 
That complete sets of eigenfunctions do exist seems plausible 
enough to justify reliance in the conclusions. 


BASIC RELATIONS 


For any one of the true solutions of the eigenvalue problem, 
for example, the eigenvalue Ca and eigenfunction %a,Pm; 
equations (41) and (42) may be written as 


r ft 1 L pL . 7 
Culo L enar L L RTP adh (A) 


and 


F TNR 
Cut lara İL Qe O ónda} (A2) 


To make the notation more concise, let the coupled mode 
Ym, Pm be represented by tm. Then if y4,Q4 is substituted 
into the right-hand sides of equations (A1) and (A2), the 
left-hand sides may be represented by Capa, Furthermore, 
because of the linear character of the equations of the prob- 
lem, substitution of the function series 


e 
>ya Wr 
m] 


into the right-hand sides of equations (Al) and (A2) gives 
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(43) 


for the left-hand sides the function series 


20101, (44) 
The coefficients a; are, in general, complex. The complex 
eigenvalues C; are assumed in the subsequent proofs, except 
where stated otherwise, to be different from each other, and 
the eigenvalue having the largest modulus is defined as Cı, 
the second largest, as Ca, and so forth, so that 


(I>e ... (Ã5) 


Expressions (A3) and (À4) are the expansions, in terms of 
the eigenfunctions and the eigenvalues, of the functions 
previously referred to as the assumed and intermediate de- 
rived modes, respectively. The subsequent proofs of con- 
vergence are based upon the fundamental relationship that 
exists between expressions (A3) and (A4). 


FUNDAMENTAL MODE 


The fundamental mode and eigenvalue are found by itera- 
tion according to the original Stodola procedure. In the 
present terminology and notation, this procedure and its 
proof of convergence are as follows: The coupled mode as- 
sumed at the beginning of the first cycle of iteration in general 
contains some component of each of the eigenfunctions; 
therefore its most general expression is 


a 
401 U — 2 a U ı 
-] 


The intermediate derived mode (which in this case is also the 
final derived mode inasmuch as no sweeping operation is 
required to obtain the first mode) is for this first cycle of 
iteration 


(A6) 


4px UU =w 2 -— Ca LU. (AT) 
fæ] 

The second and following cycles are begun with the final 

derived mode of each preceding cycle, and thus the assumed 

and derived modes of the nth cycles are 


m Ca arı (AS) 
ml] 
ap, AFD => C anpi (Ã9) 
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In accordance with the definitions given in equation (45), 
all terms on the right-hand sides of equations (48) and (A9) 
except the first are negligibly smell in comparison with the 
first for large valuesof n. In the limit the fundamental mode 
is obtained as 





lim w 0*5 =lim Cg ıt; (410) 
R> no 
and the fundamental eigenvalue is obtained from 
: wy (14-1) "n 
lim ur -vI (411) 


TRANSFORMED SECOND MODE 


The initial assumption of the transformed second mode 
in general is of the form 


-D ar ^ 
ON 
= Ur A 


in which the arbitrary coefficients 5, are in general complex 
and the subscript À refers to values of either the flexural or 
torsional components of the eigenfunctions at station À. 
More specifically, if, for example, the nodal (zero) point of 
Wa is selected to be at station À in the flexural component, 
then the subscript A refers only to the flexural components 
Of ni, Vs, Wa, - . - and not to their torsional components. 
Thus each term of the series in equation (A12) satisfies the 
requirement that either the flexural or torsional component 
of the assumed mode be zero at station A. 

To simplify the subsequent work as much as possible, the 
eigenfunctions are henceforth assumed to be normalized to 
unity at station A; thus 


(412) 


(sa = 1 ü= 1, 2. Ölə ) (A13) 
Equation (412) now takes the simpler form 
a => b,(w,—wi) (A14) 


The assumed mode given by equation (Alt) lends, ac- 
eording to equations (À3) and (À4), to the following inter- 
mediate derived mode: 

m= b ,(C Op — Cw) (415) 
Sweeping of this intermediate derived mode with the first- 


mode shape (previously determined) leads to the derived 
transformed second mode of the first cycle as follows: 





Ga = py — =>) Cow) (418) 


When each succeeding cycle is begun with the derived 
transformed second mode of its preceding cycle, the vartous 
functions for the nth cycle are 


Un" = Cb wu) (A17) 
w,™ = 23 Catb (Cw One) (A18) 
270 = Ci bao; — wy) (119) 
The limits as n approaches infinity are 

lim Way AFD = lim Û ,*bı(at'ı-- iw) (420) 

nov no 

and 
a w in T1 

ne “əə | - 1) 


Equations (420) and (421) show that convergence to the 
exact-transformed-second-mode shape t—w, and to the 
exact second eigenvalue C, can be obtained theoretically. 


TRUE SECOND MODE 


The key to computation of the true second mode is readily 
found in the simple case illustrated in figure 1. In this case 
the sweeping function of the final cycle of iteration would be 
the displacement produced by the forcing load + (w?—w,”)y,, 
in which y, is the first-mode component of the transformed 
second mode Ya. The sweeping function is designated by 
Ya Which has a well-defined numerical value in the iteration. 
Thus the value of y, could be found from the equation 





yla — os (y? 
Ua m n. (s 1) Hy (A 22) 
that is. 
— Ya 
dm (A23) 
I 
Q1 


The sum of Ya given in the iteration and y, given by equation 
(A23) gives y», the true second mode; that is, 


Ya TY Va— Vi t Vi Vs (A24) 


By analogy, the true-second-mode shape in the general 


(complex) problem under consideration is found as follows: 
The limiting value of the sweeping function is, from equation 


(418), 
— an, 02 
= 1 w=lim 5 (2-1) baw, (425) 
Cy 





lim 0 =lim 


naz 
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The expression analogous to equation (A23) is 


lim 95,” 
CO =lim Co" bw (A26) 
O, | 
The expression analogous to equation (A24) is 
lim wy, 
lim Wat += D Cy baw (427) 
N—0 Boca! SEN 1 N—+00 ; 
Cs 


which gives the exact shape of the true second mode. 
TRANSFORMED THIRD MODE 


The first cycle of iteration for the transformed third mode 
begins with an assumed mode that has two zero values, one 
of thege being in the same (flexural or torsional) component 
and at the same station (station A) as previously employed 
for the transformed second mode. The other zero value 
may be taken in the same component as was the first zero 
value and at a different station (station B), or it may be 
tuken in the other component at any station, including 
station A. Either of these possible selections for the loca- 
tion of the second zero value is indicated in the following 
equations by use of the subscript B. The initially assumed 
transformed third mode may be written as. 


as = > d, | o —. — (ou) (A28) 


1Ds— 





in which the arbitrary coefficients d, are complex. Each 
term of the series in equation (A28) is zero at station A by 
reason of the normalizations stated in equation (A13), and 
cach term is also zero at station B. 

The various displacement functions for the general (nth) 
cycle of iteration may be expressed as follows: The assumed 
mode is 





aa = Y C, 71d, | vim — (= 0 — wı) | (A29) 


— 


The intermediate derived mode is 





wy Z C," 71d, | Com C. (E 7 əl (Cap — Coo) | 
(A30) 


The result after sweeping the intermediate derived mode 
with a, first-mode shape such as to make the sum zero at 
station À is as follows: 


ap, ə) (Dem ək C, 71d İc 
w= >; t i| Cı@;--ı- 


ə) Cs) | A(81) 
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Sweeping of the mode given by equation (A31) with a 
transformed-second-mode shape such as to make the sum zero 
in the flexural or torsional component (as the case may be) at 
station B gives the derived transformed third mode as follows: 


(5) 
aD 
py 7) — 0 TA 
+1) ,__/ Metn Wi Za 
Was Mp — — ) w| —————”” ¿(002 01) 
101 Za 104—170, 
=> Cedo | oo (= m) (oe | 
w— 


The limits as n approaches infinity are 





(A32) 





+ — 
lim Wa =lim Coda | um -( : 


n— Wa— 


Bt) (o) | (A33) 


and 





(A34) 


(8) L 


As shown by equations (A33) and (A34), convergence to the ex 


pə), (w— wr) 
and to the exact third eigenvalue C; can be obtained theoreti- 
cally. 





act-transformed-third-modeshapew,—w; — 


TRUE THIRD MODE 


Computation of the true third mode is explained by refer- 
ring again to the simple problem of pure flexural vibration 
in which air forees are excluded. "The transformed third 
mode in this simple problem would be given by 





Loc. do Nu Wl X ; 
Mans Ua Y (11) Gn y) (A35) 


The total load required to hold the beam in equilibrium in the 
shape ya is 


rata? | 1 (E pt) y yu? A y, (430) 


If the beam is vibrating with shape yış at frequency ws, the 
inertia load is given by 


¥ Ws ya yoq” [ v- (Dto ya vı) | (A37) 


The forcing load required is the difference between tlıe total 
load (expression (A36)) and the inertia load (equation (A37)), 


that is, 
y (og? — w’) 6 nt) | Vi tas” 





(A38) 
The displacement produced by this forcing load is 


0... 
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and this displacement must be equal to the sum of the 
sweeping functions in the last cycle of iteration (if the 
iteration has been carried to complete convergence). The 
first sweeping function is of the first-mode shape and the 
second sweeping function is of the transformed-second-mode 
shape. If the expression (A39) is written in the form 


= p- pe Gi) (LEB) at 
— yn 


(57 


euch of the sweeping functions contained in the displacement 





(A40) 





produced by the forcing load is obvious. Thus 
(96 Ja UL X 
ə ei PRE: m) (s -1) hA y, 
(A11) 


and 


Jus ( == 1) (Py) yü (A42) 
in which yər and yas designate the first and second sweeping 
functions, respectively. Both of these functions have well- 
defined numerical values in the iteration. 

If now a simpler notation is adopted, equations (441) and 
(442) can be written as 


Yo: =( 2 — (lou E 1 nə (A13) 
atid 
Ym=( E 1) Je (A 44) 
in which 
n= 1225) ME o 
und 
ni” a Tem y) (447) 


The true third mode is clearly given by the sum of equations 
(435), (445), (446), and (A47); thus 


Mja Yu Vis T Ves (A48) 


The transformed third mode Ya is given directly in the 
iteration. The procedure for finding the other components 
ou the right-hand side of equation (A48) is as follows: 
Component Ya, by equation (444), is 





(A49) 





212483—54--.—22 


Component ji; is known when Ya is known because its rela- 
tion to Ya was established previously in connection with the 
transformed-second -mode calculations (see equation (A24)). 
Component zr is then found by equation (443) as 


(A50) 





By analogy with the foregoing case, the true third mode in 
the complex-eigenvalue problem is found as follows: The 
limiting value of the second sweeping function is (see eque- 
tions (431) and (A32)) 


(x) 
UN 
P0,''-- 5 ) mi 
i e! ---- lim] —— 2 | (m—uw:) 
no Ra Wai— Wy, 


lim aaa 
—lim €,” (&-1)4 an UA) (105— Wy) 


The limiting value of the first sweeping function is (see 
equations (430) and (431)) 


10 (a) 
b ~a) ib 
n-—u 


lim wa — — lim 
ENT š Cı Cı Cy Wi— WD: 
“lim Ced, LS LLE TLE Lle 


(A52) 








(451) 








The quantities analogous to yı, and Ya of equations (443) 
and (444) are, for the present case, lim w,.‘* and Itm wy». 


nov 11-+ co 
The latter quantity is obtained from the relation analogous 
to equation (449) as follows: 





Lim Was" —ap 
lim 17, =" iia Ced (— t) (mw) (A53) 
to 2. 1 To 
C. 


The relationship of lim 1,2 and lim w,,“' is obtained from 
n- m na 


equations (A26) and (420) of the section dealing with the 
transformed second mode. Thus, 


lim uy 
lim 77,2” =lim e Jal Er tasa | azə 


Faren ne [(Wa t * Ps] aş. (420) | Cro ¡ 
2 Eg. (425) 








(A54) 


The quantity analogous to zu of equation (443) is, for the 
present case, lim :3, and is obtained by an equation 
no 
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analogous to equation (A50) as follows: 


C, 
lim Vw — o i lim wu 
8 , 


lim Uy U) === 


HR— co 





=lim C,*d, (a — 


no 





The exact shape of the true third mode ws às given by the 
sum of equations (A33), (A53), (A55), and (A54), which is 


lim (tw; "t? + Waa +01 tus =< hrm C "dyap a (A56) 


FOURTH AND HIGHER MODES 


Extensions of the proofs to modes higher than the third 
can be made in a manner similar to the foregoing proofs. By 
this means, the iterative transformation procedure can be 
proved, under the assumptions stated at the beginning of this 
appendix, to be convergent for all modes and eigenvalues. 


CASES OF EIGENVALUES HAVING EQUAL OR NEARLY EQUAL MODULI 


For a representative case, suppose that 


lCıl>iCal; 1611 04»..: (A57) 
and that 
[Ca — [Cs (A58) 
or that 
[Cal = |C] (A59) 


Under conditions (A57) and either (458) or (459), the 
assumed and derived modes after a few cycles of iteration 
will be virtually as follows (see equations (417) and (Å19)): 


Q 71b, (ap | — w1) + C," 5, (w — w) (A 60) 
Wa 7 * D = C3" bs (We— un) + Cı*bı (apa — w) (A61) 


Weg ™ es 


If |C,| is only slightly greater than [€;|, the second terms on 
the right-hand sides of equations (A60) and (A61) become 
negligibly small very slowly as n increases, even though 
they do become negligibly small as n approaches infinity. 
If [C;| and |C;| are equal, these terms never become negligi- 
bly small. Thus, the problem of circumventing this slaw 
convergence or apparent lack of convergence arises. 

A satisfactory method for coping with these conditions is 
to combine linearly the results of the last two cycles of the 
series of iteration cycles that have been performed. For best 
results in an actual problem, not less than the third and 
fourth cycles should be used for this purpose in order to 
reduce as much as practicable the effects of all higher-order 
components. 

The following formulas for combining the results of the 
last two cycles are based on the assumption that the assumed 


and derived modes in each of the cycles contain only com- 
ponents of the types in equations (A60) and (A61). 

The two components (with shapes w,—w, and w,—u,) 
clearly appear in the last cycle in proportions different than 
in the preceding cycle. (The proportion in each cycle is a 
complex function of the spanwise coordinate.) Because of 
this differing proportionality the results of cycles n—1 and n 
can be linearly combined so that the combined functions 
contain only one of the components w,—1, and wy—w). 
Accordingly, the ratios of both the flexural and torsional 
components of the combined functions at all stations should 
be equal to each other. In algebraic terms, this statement 
means that 

(mm Es), =R (A62) 


rw (n —1) +wv4 (1) 


in which 7 and R are (complex) constants, and the subscript 
S designates that the ratio may be evaluated at any station 
S, that is, that R has the same value for all stations. AH w 
functions must be the same type of component, cither flex- 
ural or torsional. 

Since S can be any station, the equality 


papa) + wy tD rap a E, TH 
Es D Hwa a) = (Çele nl “Duas =), wee) 
exists, in which stations 1 and 2 must be different or may 
be the same, depending on whether the w functions on the 
left-hand side are the same or different types of components 
than those on the right-hand side. The two values of r 
that satisfy equation (A63) are 




















4D Eb ATT RTS T Acar 
po ky Sena )— gema (464) 
in which . 
(Wa t) (ga, 
A-D.) (465) 
(Was), (Waz ™)e 
(Wa), (Was 95 
APAGA — (A66) 
(watt?) (Was 0*9), 
(Wa 0D), (Das 79), 
(2406-1054 — j i (A 67) 
(Wwa tt”) (pa 0*0), 
The corresponding values of R are 
ASTD n+ A (0-», At) 2 A or D 
CO ¿ED DAPA / A46-Da 
(n), (x1) (n—D.(m) 
2007007 (A68) 


r 
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These values of R are equal to C, and Ca and the corre- 
sponding values of r, when placed in the expression 


Pas (x) Wo TD (A 6 9) 
give modes of the shapes u,—w, and uy—w,. When [Ci] is 
nearly equal to |Cj|, the appropriate set of R and r to give the 
lower transformed mode w,—u is evident. When |C%| and 
|C;| are equal, the mode obtained by equation (469) with 
either value of z may be used as the transformed second 
mode, but the trends of the eigenvalues that have been or will 
be determined at other values of the reduced frequency k 
may be used as a guide in making the selection that fits the 
trend. 

In actual computations, one further cycle of iteration be- 
ginning with an assumed mode given by expression (A69) 
should be carried out to assess the extent to which the fune- 
tions 45779, wee, and wee*t® are free of all except the 
two components of the types appearing in equations (460) 
and (461). If the ratios of this cycle are not reasonably 
constant, the unwanted components still present have to be 
removed by carrying out another cycle of iteration and 
again applying equations (A64) and (A68). 

The method just described is clearly applicable in the 
general cases |C,{=[Ca4:| or [Cal = Cal. 

Eigenvalues having equal moduli include the special case 
of identical eigenvalues. As a basis for discussion let it be 
assumed that 


¿Cu> Cal zi c6... (A70) 


and that 

Cass C, — Ca (A71) 
The significanee of the occurrence of these two identical 
eigenvalues is that the wing system may oscillate with the 
same frequency and artificial damping in any of an infinite 
number of modes, any two of which are linearly independent 
of each other and of the first, fourth, &nd higher modes. 
This infinite number of possible modes (all corresponding to 
(4; are the infinitely many linear combinations of two 
basic linearly independent modes that are necessary and 


sufficient in combination with the first, fourth, and higher 
modes to describe an arbitrary displacement of the wing 
system. Clearly, only two linearly independent modes cor- 
responding to the double eigenvalue Ca are required for 
analytical purposes. These two are designated w, and wy; as 
before but with the reservation that us and w, must be 
derivable as two differing linear combinations of a single 
basic pair of linearly independent modes that also correspond 
to Cos. 

Equations (A20) and (A21) are replaced in the present case 
by 
lim 10,,*? = lim CA” 152010: ri) + baw — wp] (A72) 


11— OB 


and 





RED 
lim Yao, (AT3) 


1-+ co Waz 


Equation (A27) is replaced by 


> lim eu ™ 
lim qp a tl E NUE s lim Co" (bw. + baba) (A74) 
"cow — EM 1 uo 
Ca 


The transformed second mode (equation (A72)) is in this 
case & linear combination of the first three eigenfunctions, 
and the so-called true second mode is actually a linear 
combination of the second and third eigenfunctions. 

If the iterative transformation procedure is now applied 
in the regular way to determine the transformed third mode, 
the third eigenvalue, and the true third mode, the results 
will be as follows: The transformed third mode will be, like 
the transformed second mode, a linear combination of the 
first three eigenfunctions but will be linearly mdependent 
of the transformed second mode. The so-called true third 
mode will be, like the so-called true second mode, & linear 
combination of the second and third eigenfunctions and will 
be linearly independent of the so-called true second mode. 
The results will also include a second determination of the 
double eigenvalue Cy. It may therefore be concluded that 
the iterative transformation procedure is valid and sufficient 
in all cases of eigenvalue multiplicity. 


APPENDIX B 


THE COMPLEX STIFFNESS FOR BEAMS WITH STRUCTURAL DAMPING 


The familiar concept of a complex force K(1+-2g)8 in 
simple (one-degree-of-freedom) vibrating systems having 
structural damping may be easily extended to continuous 
vibrating systems such as beams and airplane wings. The 


quantity K is the elastic-spring constant, s is the displace- 


ment, As is the elastic-spring force, and Kgs is the structural- 
damping force. 

For a beam in flexure, the stiffness of ihe fibers is given by 
the modulus of elasticity E, which is analogous to the quantity 
K for the spring. The elastic stress at any point of the cross 
section is given by «E where eis the strain which is analogous 
to the displacement s. Then the complex stress at any 


point of the cross section of a beam with structural damping 
is E(1+ig)e. The complex bending moment corresponding 
to this stress, obtained in the usual way by integration of 
the moment of the stresses over the section, is EI da lar zı ə 
This result leads to the concept of a complex — 
ET(--ig) for beams in flexural vibration with structural 
damping. Similarly, the complex stiffness of beams in tor- 
sional vibration with structural damping is GJ(1--9,). 
The subscripts y and ¢ indicate that the structural-damping 
coefficient g may have a different value for torsional vibra- 
tions than it has for flexural vibrations. Both g, and g, may 
be functions of the spanwise position z. 


APPENDIX € 


FORMULAS FOR EQUIVALENT CONCENTRATIONS AND INCREMENTS OF TORQUE 


The formules used in the numerical examples for com- 
puting equivalent concentrated loads and curvatures are 
those that have been derived in references 7 and 8. For the 
concentration at an end station the formula is 


a (7p, + 6p.— Ds) (C1) 
At an intermediate station 
p: = ig (pit 10p.+p3 (C2) 


The significance of the quantities used in formulas (C1) and 
(C2) is shown in sketch 1. 


Do - -— 
= 
-À 


nunun 


A 


Bketch 1. 


Distributed -Joad curve 





These formulas are based on the assumption that the 
distributed-load (or curvature) curve is a series of second- 
degree parabolic ares. When applied to distributed flexural 


loads, the formulas give concentrations which produce the 


same bending moments in the wing &t all the selected stations 

us the distributed load. The formulas may be correctly applied 

to distributed torsional loads only if GJ is constant over each 

bay. In this case the formulas give concentrations which 
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produce the same torsional displacement at all the selected 
stations as the distributed load. For a station placed at a 
discontinuity in ordinate or slope, formula (C1) must be 
applied to both the left and the right of the station and the 
results added. 

The formulas for obtaining increments of arca beneath a . 
curve of distributed torques are derived in reference 8. 
These formulas are based as before on approximating second- 
degree parabolas. They are given here in a slightly different 
form which is better adapted to present uses. Thus 


Aıê (9144927 04-5 (i—qə (Ca) 


A=) (qi--44:--g3—? (91— 13 (04) 


where the significance of A, and A, and of q, ge, and gy is 
shown by sketch 2. 


increments ofc... 
` _ 
torque 1 š ^ 
E j \ 


_- Distributed- foad 
di curve 





Sketch 2. 


The ordinate at a discontinuity should not be used as the 
middle one of the three ordinates selected for use in formulas 
(C3) and (C4). The formulas are valid only where the three 
ordinates are connected by a continuous curve. 
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TABLE 1.—ITERATION TO OBTAIN FIRST COUPLED MODE 


[Common factors for each column are 


Z | 
9 iS EKE KAKA - <ar ea us —^" 
Fd 
5 4 3 2 l 
Station 


Flexure: (a) First cycle 






aW TÛ, A my.” 


6 è 


Station 






0 
14, 30 






: = ı7. 95 | 0 | 2451 47. 86 
| | 58. 07 31. 82 
1. 13 1. 13 79. 68 
| 59.20 | 32.45 
112. 13 












Station 
Ag y ol 
Ek 
1 0 1. 397 0 1.397 | 0.594 | 0.594 ê 0. 59 
2 0 . 794 0 . 794 . 800 . 800 l . 80 
3 0 271 0 271 -219 | [219 | ıa 68 0 —14 41 
4 0 . 076 0 . 076 . 086 . 047 





AN ITERATIVE TRANSFORMATION PROCEDURE FOR NUMERICAL SOLUTION OF FLUTTER PROBLEMS 


FOR be FOR WING WITH CONCENTRATED MASS 


given under the column headings ] 


ET — — .—-9 
Zİ 
5 4 3 2 | 
Station 
(a First cycle— Continued (b) Second cycle 


| i 

























(c) Third cycle 
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eng | t i Í 
14 1:15 ! 16 | | IS ` 19 3 : 1 2 
| | | | ! | | : 
| | | | | : 
" : = y yı") yi Lyn 
a | a | B | Ay yA | yid 3 y; no ya yi ® Xy 
00 en | : 
MA by hathy Mtby | Mi hiy | datby Roly 
"Ple S | ph" 0” |. ER" Ela” | P 1 EI," 
280. 27 | 280. 27 281. 65 1 000 | 281. 90 | 281. 90 
120. 90 | 120. 90 
15.90 j| 15. 90 159.37 | 280 232 . 570 | 160. 36 | 282 
: 108. 00 | 108. 00 
(17.53 İf17. 53 
71.46 | 39.15 
79.75 | 43 70 15. 22 | 282 284 . 055 15. 33 | 279 
27. 76 | 15.22 
50. 60 | 27. 76 0 0 0 
(a) First cycle—— Continued (b) Second cycle (c) Third cycle 
m 
dı) | yer 
$0 La? 
Ae 
EIL 
, —12 28 | —1 662} .... — 0. 00594 282 
„ 5 
—12.87 | —1 742} .... —. 00622 —1. 861 299 282 
1. 
— 14.26 | — 1. 930 5 —. 00689 — 2. 048 297 282 282 
—7. 15 TAL! —.963 | .... —. 00344 — L 022 297 282 
—7.11 0 0 0 








EI, 
My 423,000 
m= V sə” “82 =38.7 radians per second 





St tion 


TABLE 1,—ITERATION TO OBTAIN FIRST COUPLED MODE FOR k= «o FOR WING WITH CONCENTRATED MASS—Concluded 
Torsion (with GJ variable): (a) First cycle 


nn —-——-———-——-.—.-.-.--.————————..X.x.X a a i P I Y Ur — rr nr 





Mü" ə l 
A 


0 0 





1,089 | 1. 080 


1.080 | 0 
15-ə do | t 
nia | qu. H 014 —14, 88 0 Uu 08 
, 104 . 090 | 


1. 704 — 14. 68 
| 1, 719 — 14. 68 


. 028 „ 015 


a 
I 
| 
! 
I 












SGE 


SIILAVNOHAY TOZ AMLLINAOS AHOSIAGV 'IYNOLLYN—8401 LHOdAHH 


TABLE 2—ITERATION TO OBTAIN TRANSFORMED SECOND MODE FOR k= e= FOR WING WITH CONCENTRATED MASS 


[Common factors for each column are given under the column headings] 


” 


2 e E x. = ves O, AA 

4. 

4 

b 4 a 2 1 
Station 


Flexure: (8) First cycle (b) Second cycle (e) Third cycle 


— — 
—. nn nya ee — 





— — oe n 





em e w -e o 









Station Ya» y, yu? y ax? 
dolby Müby 
ER” b . ə 
—411.0 | 411, 0 0 —465,8 | 465.8 | 0 
— 252, b | 234. 2 3,8 İ 26,4 —.820 | —285. 4 | 265.4 | —20. 0 
—90, 7 | 80.0 | — 882 | -1ı2: 0 | 90.0 | —21, 4 
—80,2 | 224 ê. û — gap | —3328| 2h4| —&4 
0 0 0 0 0 















Station 
2 — 1.15 | 12, 28 , D05 26, 24 i ; 
d — 1, 28 | 12, 12 „ 988 26, 40 | —3, 40 | 23, 00 | 23.3 23, 8 
4 —, (13 0, 07 408 ıa, 20 | —1. YO | 11, 50 | 23,8 
0 0 0 0 0 0 0 














e RE 27 wî; w= 188.8 radians por second 
Aa naz A 


SKSISONd WALLATA 40 NOLLNTOS TVOIHHNAON HOZ WHügWoOHd NOILYRHOISNYIL WALLVHWILI NY 
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— — = sus -ٌ 


TABLE 3.—ITERATION TO OBTAIN TRANSFORMED THIRD MODE FOR k= < FOR WING WITH CONCENTRATED MASS 
{Common factors for each column are given under each column heading] 


Neat ee 85... 













































2 ; 
b 4 3 2 1 
Station 
Flexure: (a) First cycle (b) Second cycle (c) Third cycle 
1 4 5 6 7 
O | Xya* 
(D (0 (9 a (3) co | Ys Ya 
Station Ya Yeas Yas Yas Vu 1 yasaq 9) | ya? 2.8 | Eya 
| b 
| : 
448. 11 — 448,1 0 0 
1. 000! 242. 6. O 
1. 000] 104. 01 — 87. 0l -- 11. 9 . 903) 91. 7.2 7. 07 
2 31.0 —24 4| .--4, 7 . 883| 20.9 7. 1 
0 0 0 0 
Torsion: (a) ‘First oyole | (b) Second cycle | (o) Third oyole 
1 2 3 4 5 6 I 2 3 1 2 3 4 5 
da? | 
da? | ds 16930) | a | der’ du da? | ət) | $490 $a | bo | $49 | evnig 
Mty Mty 
EI, Br o 
— 15. 79| 2, 84| 12. 95 0 0: —9. 4 2. 57 6. 89 0 0 — 9, 46) 2. 56 6. 90 0 
—. 083| — 15. 81| 2, 96! 12, 88 . 13 | —1,8 02521 —9. 50) 2,08 6. 85 03 L 3 . 0042] — 9. 50; 2. 68 6. 86 . 04 10,0 
—. 3041 — 18. 03] 3. 26] 12. 80 . 03 —, 1 00541 —9. 721 2,96 6. 80 04 7 „ 0056] —9, 72] 2.94 6. 81 . 03 5 7.15 
— R, 001 1.63] 6.40 . 03 —. 2 0054] — 4, 85) 1.48 3. 40 03 6 . 0042] -- 4. 851 1 3. 41 . 03 7 
0 0 0 0 0 0 0 0 0 0 0 0 0 





l(XEya?,X4a9Y 711 MÜ. ce a= 
2 (Eye, Esa 7.11 eyy Pi 01244 radians per second 


OE 


SDILNVKOHHY HOZ HELLIANOS XHOSLAQY 'IVNOLLYN—@201 LHOdAHH 


TABLE 4.——AERODYNAMIC-INERTIA FORCE COEFFICIENTS FOR VARIOUS VALUES OF k FOR EXAMPLE WING 


[Common factors for each column are given under the column headings] 


Flexure 





SIXWISONd HELILNTA JO NOLLOIOS IFOJHARAN HOS WHulquoO0Hud NOLLVTYUOQASNYEL AALLVURLI NV 


TEE 
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tig 


Station 


. 028 
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TABLE 5.—ITERATION TO OBTAIN FIRST MODE FOR 


[Common factors for each column are 








£ 
4 
5 4 3 2 1 
Station 
Flexure: (a} First cycle o 
2 8 4 5 6 
tir) Pasir Prepr Prsî/ı Pied 
y 2 
H mem 
0 3L 0 0 0 0 
0 13. 77 0 
0 3. 44 0 0 "n 
0 868 0 D 
0 0 0 0 
Torsion: (a) First cycle 
1 2 3 eub. s 5 6 7 8 
ır ® du) Qrsyr Ürer Qr,yr Uad: da qr 
bêv m A,biy q? 
H A 
0 0 2.30. 0 0 2. 36 0. 039 
Ü 1. 048 0 0 0 L 048 1. 092 | 
0 0 282 0 0 0 . 262 { vr 
0 0 . 0660 0 0 . 0660 ` 077 
0 0 0 0 0 0 0 
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4£=0,1443 FOR WING WITH CONCENTRATED MASS 
given under the column headings] 


Station 
ta) First evcle—Continued 








iran — n. — -— r X— — ——. ——. 
` 


I7 18 19 20 
e a B Ay 
Marby 1 Agby I hotbry 
Elu C Ela € Ela 
95.3 95.3 
13. 5 13. 5 
8L 8 81 8 
14. 5 ( 14.5 
{ 23. 2 12. 7 
54 6 30. 0 
6L 5 33. 7 
20. 9 -IL5 
81 03 38. 2 20. 9 
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TABLE 5.—ITERATION TO OBTAIN FIRST MODE FOR 
Flexure: (a) First cycle— Continued 












23 2 2 31 
r En m 
Station Pau Ter —Prnyn Pr PO, 
— 4. 31 —4. 31 
— 6. 00 — § 00 
— 1. 099 
—.-234 0 0 -—1, 333 
—; 198 —. 103 


Torsion: (a) First eycle—Continued 





18 19 20 — 21 29 23 şî je 
Station A Özer — Qıslin | —Qrér qı qı q Ovi 
Ey w | f | MET et ably wot 
" j A A 
4. 05 1. 612 1. 612 
1. 798 1. 874 1. 874 
. 412 . 412 
0 . 450 { 0 - 412 ü 
0 . 1143 . 133 . 073 
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k=0.1443 FOR WING WITH CONCENTRATED MASS—Continued 
(a) First eycle—Coneluded 












| 32 | 33 | 34 | 35 | 36 37 | 38 | 39 | 40 | 41 | 
| 0 7 00 e DEE — 
” . = | a | co gua? + iyu? 
| Vr | AM; | Mr; | a | o | a B | Ay | Yu | yix 9 + iy; P 
i 
My | My a 
ü T : | i 
T¬¬ | | | 
—4 31 —4. 81 | 
| —4.31 | | | 
| —9. 31 —9. 31 | 
— 18. 62 | i i 
| - 10.64 | —5. 84 | | | | 
i | —19.46 | —19.46 —19.46 | —10.67 | | [| sen BQ culos 
` —1083 | —85 98 | | | —6.41 —3.52 | 
| —25. 41 | —25.41 ,  —11 70 — 6. 41 | | | 0 | 
t , ! 
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Station 





Station 


yır? 







— 0. 00382 
—. 00407 
—. 00457 
—. 00228 


É 
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TABLE 5.—ITERATION TO OBTAIN FIRST MODE FOR 


Flexure: (b) Second cycle (c) Third cycle 
5 2s 1 2 8 4 
PNA l 
ELM yis yii yia? puto 
iy 1 


a ei yum 





(e) Third cycle 


1 2 3 4 





bed “idar” 





(3) (s; (0) (4) 
¿e hip dz f gi də di 

Mir 1 

Elp C 
0. 00466 : i j .| --0. 00680 | 0. 00399 — 1. 570 1. 594 
. 00358 000" İ nit . —. 00680 . 00300 — 1. 690 1. 350 
. 00126 Š —. 00716 . 00068 — 1. 963 . 764 
. 00064 . 912 2002. | a . —. 00358 . 00036 —. 082 . 887 

0 00100 : 0. x 31-70 0 0 


LF Z (vına + mır) Z (oia? + iÓ rh) 





E XY 1. /423000 ; i 
3 Y Gg in y a |= (269.5 — 82.26) Ela C w == “360.5 7799.7 radians per second; 
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k=0.1443 FOR WING WITH CONCENTRATED MASS—Coneluded 


te) Third eycle— Concluded 








CT: 
| | 
tre + tyr (9 o 
| pir + iyi? La a 
Mty 1 | b 
Ela C 
267.5 — 82.890 | L 000 0 
25 | - 569 . 0030 
268 — 83.21 | . 194 . 0027 
e se aê i - 054 . 0011 
| 0 


(d) Fourth cycle 





3 4 5 6 ` 
yın? gu yis + rr E (ns + iyu) 
yia? + 1yu i na + ty) 
Xby 1 Aşiy 1 
Ely C Ela C 
268. 6 — 82. 56 268.6 — 82.561 
153. 0 — 46.29 fk baa 
52. 4 — 15. 28 | 270-- 82.5: 269 — 82.81 
E^ ré E iN. do eee owes 
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fred + is (4) (p a o dia + ig? i. (e: O + ipur) 
| dig? idi) | a i in is dig? + (du L (pir + ipu”) 
| Ay I 
BT, C 
| — 0. 00704 0. 00379 — L 591 L 589 270 — S0.8: 

| —. 00720 . 00281 — L 711 L 340 | 5 
282 — 79.81 —. 00751 . 00052 —L 977 „ 758 269 — 82.47 270 — SI.5: 
MERMA | | — 00376 , 00028 — 988 1384 j| .......... 
_ —8$22_ „n 1 39.7 _ 

ga 289.5 =— 0.305; = 3 014437916 feet per second 
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TABLE 6.—ITERATION TO OBTAIN TRANSFORMED SECOND 


{Common factors for each column are 


Z 
2 | e | l e B 
5 4 3 2 l 

Station 


Flexure: (a) First cycle 


Var? + iy aÑ 
Vos D + eyar) 


Station 


m m Amo m m eee ee = s m, = A xü = 





1 2 3 
- bar? + idu? 
(0 o o Par” Tider 
Station Tat du du 7 bar + 19 4,0 





en Hn Q2 O ked 
cn 
bo 
emo 
0400 
— 
uu 
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MODE FOR k=0.1443 FOR WING WITH CONCENTRATED MASS 


given under the column headings] 


A ! e ! | 

7 0 SERE es at 

Z 

5 4 3 2 I 

Station 
(b) Second cycle 
ST. ft O n | | 
1 2 3 | 4 5 6 7 & 












D + £y ua; 

"T | AD e ra ra (n a a Zen “Zer, 

UatR Vaezi HER | Yar | VitR | Yar Yer Year Yar + ives: 
k 






| b 

0 | o 

—. 188 ' - 064 

Sle 3 . 090 

| < Oi |! . 036 
| 0 

(b) Second eycle 
l | 2 | 3 4 5 





Parr + ipa 
bær? + bar? 


| 
| 






- — 


l 
I 

—. 0426 | 21. 60 — L. 59 —6. 61 
| 


1. 000 30.80 | -e92 | -5 105 24. 92 —2. 87 21.92—2.87i 
` 886 —.0153 28. 36 —& 33 —& 12 3 21 22. 24 —312 ea 
E i | 124 14 99 -23 Í 262 zs | 
` 298 — 0221 İL 01 —2 36 —330 | 163 TU 557 15 

0 0 ıı | 0 | | 
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TABLE 6—ITERATION TO OBTAIN TRANSFORMED SECOND MODE 


Flexure: (c) Third cycle 









1 2 3 4 5 E 6 8 9 
A | —— € Rə S 
2+ (0 L fy 4 (0 
(3) (b a) (8) (0 (5 (2 Van Yas! 
Station York Yazi Yer Ç Yar Ung D Nat Ya F dua 
| Mby 1 wri | 
Elp C BI u O 
0 0 —931.9 | 156. 25 931.9 | 215625) 0 0 | 
—,271 . 149 — 537. 2 92, 18 530. 5 — 86. 3 — 6.7 5. 9 28 — 6.31 
—. 190 . 161 — 187. 2 32, 98 181 2 — 27.7 — 6. 0 5. 3 31— 1,4: 
—.071 . 065 — 58. 0 9. 60 50. 8 — 7. 41I — 2,2 2. 19 32— 1.41 
0 0 . 0 0 0 0 0 0 
Torsion: (e) Third eycle 
| 1 2 4 ; 
. DL O 
s bare da PS par? + bari" 








bart + thar 









20.38 — 3.82 









29. 90 —7. 81 — 6. 29 3. 71 26.3— 3.00: 
609 —. 0644 25. 10 — &. 99 — 6. 03 1 64 27.0— 4.29 
314 —. 0333 AL 80 -à 08 —3. 46 » 84 27.0— 4.26; 


daly 


E (arr ® + idar ®) 





1 [ Z (yar + tare) |= - 3. 
2 LE (Yare + ty cer) + E (bar Hip] 7 HABEMUS een Flu C* 
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FOR k=0.1443 FOR WING WITH CONCENTRATED MASS—Concluded 


(d) Fourth eycle 











| 1 5 6 
| Yarn? Yarr? Var? 

‘ b Aotby I 

Elp € 

0 ' 0 —941.0 | 165.58 | go | - 165. 58 
| —.981 .184 | —543.3 98. 06 534 9 —91. 8 
| —.252 [165 | —190.0 35. 82 183. 1 — 29. 50 
| S : 





w= 121.0 radians per second; ge — 0.1308; r, =280 feet. per second 


Va*R 


Past 


Ver + tU az? 


Værd + Eyaz? 


Gaz” + idas 


bear Hipan” 


27.13— 4.07: 
27.1— +.05 
27.7 — 4.02: 
27.5— 41.066 


E (Yan Ligar) 
(Yer? + iyeu ®) 


31—3.41: 


L (par ® + idear O) 
MECA iper?) 


27.3— 4,051 
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TABLE 7.—ITERATION TO OBTAIN TRANSFORMED THIRD 


(Common factors for each column are 


342 


2 | | | 

Z“Zl]———u@ 

/ 

Edi 

5 4 3 2 1 
Station 


Flexure: (a) First cycle 







1 11 | 
J (1) E Yar + Pyas” 2: 
Vaan Van” +t Yaar) | 


Station 






MM mm eae meme x 





"Par + ipes” 
bare + tó tD 





Station 
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MODE FOR k=0.1443 FOR WING WITH CONCENTRATED MASS 


given under the column headings] 





a | a | 


5 4 3 
Station 


(b) Second cycle 


bə 
aa 











“m. | | 
Í 2 Dos | ı Í 3; | 6 ` = 
MINES TIMENS. o e VENE re q 
| | M4; 
o ce a | o | e | a Vosa + iya 
Hai | Ves | Yar | yər | VER E | Viazn Ysaur™ Ya Tiya +i 





—. 33 17.5—0.13: | 


0 

16. 9 

15. 3 .21 | 15.84-0.21i 
P x 









hesr ® F 14:43:78 


| 











| Paro + ipar” 
—15.90| 1640 312 | —2392 | 1287 | 0.752 0 0 | 
—16.02| 1.521 3.28 | -L988 | 1161 .358 | —L15 | —.059 | ------------- 
—16.00 | .858 ; 3.50 | —.872 8.23 | — 460 , -418 | —.474 | 17.4—0.116/ 
—8.07 | .472 | 180 —.445 123 | —.242 | —204 | —. 215 18. 0— 0.0931 
0 0 m. | 0 0 | o | 9 | 
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TABLE 7.—ITERATION TO OBTAIN TRANSFORMED THIRD MODE 
Flexure: (c) Third cycle 





Yan + £y 0 


Station Vanil +1Y0 V 


515.9 | —102.86 | —515.9 | 102 86 

311. 8 — 59. 83 | —293.8 57.3 — . . 12.2+ 0.167 

118. 5 — 20. 69 — 100. 5 18. 51 . : 14.14 0.107 
i 2 E 93 Eo 2 E 00 : s .5 | 12.7 — 0.067 





Par TL ra, rtt 
par + ipay Y 


Station 





1 [ Z Yyar” + tar) o | Xy 1 
3 LE (yas tyes) Y (Gar 9 Fo | (1482+ 0.4457) pr 5 
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FOR k=0.1443 FOR WING WITH CONCENTRATED MASS—Concluded 
(d) Fourth cycle 





| 

; ! = 

| | E 

j a ! (4) var? + f yar? E (Yard + LU ar) 

| Yate | Yan | | DT S (yea F iya] Fiyat M 

x E 
b 

| 0 0 494. 8 |-86. 24 |-404.8 | 88.24 | O 0 0 0 

| 952 | — 033 | 300.3 |—50. 13 |-281 6 | 47.75 | —5.1 | 235 | 13.6 | —.03 14. 3+0.46i 

1.000 | 0 115. 2 |—17. 33 | —96.2İ 15.31 | —43 | 253 | 147 . 51 14. 740.51 14. 5340.5471 

| .390 | — 005) 3441 —487| —26.0İ 4.12. —L7 | LO 5. 8 . 25 14 840.835 . 

: 0 0 0 0 0 0 0 0 0 0 | 

I 
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«3 = 168.9 radians per second; g. = 0.030; n= 390 feet per second 
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TABLE 8.—COMPUTATION OF TRUE SECOND MODE FOR k=0.1443 FOR WING WITH CONCENTRATED MASS 


| Flexural functions aro in terms of b: torsional functions are in radians; gi = Fi 8.05— 1.600; | 
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TABLE 9, —COMPUTATION OF TRUE THIRD MODE FOR k==0,1443 FOR WING WITH CONCENTRATED MASS 


[Flexural funetions are in terms of b; torsional functions are in radians: an = Fyr 16.97 — 6.00; GE 1 Pa < 0,940 0.3197 
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